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Abstract. We consider the general bilinear control systems in three-dimensional Euclidean space which satisfy the 
LARC and have an open control set U ; and give sufficient controllability conditions for corresponding projected 
systems on two-dimensional sphere. 



Introduction and Statement of Results 
Bilinear control dynamical systems in a Euclidean space x e i?" are given by linear equations 

(1) x{t)^A{u)x{t), (S) 

with the right-hand linear operator A{u) depending linearly on the control parameters u — (u^, u'^) G [/ C 

(2) A{u)^ A + u^Bi + ... + u'^Bd, 

where A and Bfc, fc = 1, d are some constant linear operators. Then x{t) is the solution (or trajectory) of S 
if 

(3) x{t) ^ A{u{t))x{t) 

for some control u(t); i.e., the function u{t) : [0, T] U . The positive orbit 0^{x) of the point x is the union 
of all trajectories x{t) issuing from it; i.e., such that a;(0) — x. The bilinear system S is called controllable if 
for every non-zero point x its orbit coincide with the whole space without its origin zero: 0'^{x) = i?"\0, which 
means that for two arbitrary non-zero points a;, y there always exists some trajectory x{t) of S which connect 
them: x = x{Q), y = x{T). We may "project" the bilinear system to the control dynamical system defined on 
the sphere S^^^ of unit vectors in i?" by considering an evolution of unit directions q[t) — x{t)/\\x{t)\\ of the 
non-zero solutions x(t) of E: for the unit vector q we define the vector V{q^u) to be equal to the tangent to 
the sphere component of the vector A(u)q; i.e., V{q, u) — A{u)q — {A{u)q, q)q. Then 

(4) 

correctly defines the control dynamical system on the unit sphere which we call the projection TP"^ of the 
system E. This definition is justified since for an arbitrary (non-zero) trajectory x{t) with the control u{t) 
its central projection x{t) — > a;(t)/||x(t)|| to the unit sphere is the trajectory q{t) with the same control of 
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the projected system E^"'. If in addition some linear operator A('u*) has all eigenvalues with negative real 
parts, while some other A{u*) has all eigenvalues with positive real parts (i.e., the first one takes eventually 
an arbitrary big sphere inside the unit sphere, while the second - an arbitrary small sphere outside the unit 
sphere), then the controllability of E is equivalent to the controllability of its projection E*"": in [CK] (see 
Corollary 12.2.6) the authors the claims that the system E is controllable if and only if 

CKl) their Lyapunov spectrum contains zero in its interior, and 

CK2) the projected system EP'" of E is controllable. 

Note, that from the linearity of E it follows that E^"" is invariant under the central symmetry q —q, 
and therefore, we may further reduce it to the system on the projective space RP'^~^ = S"~^/Z2. Wc do 
not do this, since for some symmetric and non-controllable systems on the sphere their further reductions to 
the systems on projective space are controllable; i.e., this further reduction does not preserves simultaneous 
controllability. It is, also, more convenient for us to work on the unit sphere in the Euclidean space. 

Next question is: when the projected system is controllable? The two-dimensional case n = 2 was considered 
in details in [AS] and [BGRS] where the authors show that the projected system on the circle RP^ satisfying 
the LARC is controllable if and only if for some parameter uq the matrix A + uqB has complex eigenvalues.^ 
Remind, that the LARC - the " Lie Algebra Rank Condition" is satisfied for the control dynamical system E at 
the point x if the evaluation at this point of the algebra Lie £A{T,) generated by all vector fields A{u)x coincide 
with the whole i?". Easy to see that the systems E and E^"" satisfy or do not satisfy LARC simultaneously, 
and that for n = 3 the system E^"" on the unit two-dimensional sphere 5^ C i?^ satisfies LARC at q if and only 
if there exists two linearly independent vectors V{q,u) and V{q,v) at this point (our CCl condition below). 

In general, the controllability of the bilinear system (1) corresponds to the transitivity of the Lie algebra 
generated by the right-hands A{u). Classification of such transitive actions is obtained in [B] (see also [BW], 
where the algorithm is given for the practical verification of this property). Our approach is quite different: 
instead of checking the transitivity we study directly the chronological products of the control system (1). The 
reason is the following. In [Ml] we considered quantum systems with Hamiltonian depending on a parameter 
u. Using the Dirac mechanism of quantization, it was suggested some new description of the Berry-phase 
geometry for such systems. The problem of measurement of such systems naturally leads to the investigation 
of dynamical systems with the control u. These systems are defined on "level" sub-manifolds M of the 
Hamiltonian in the space of parameters; and are not linear systems. We verify the controllability of such 
systems with the help of the "closed orbit controllability criterium", sec [M2], which works for arbitrary non- 
linear control systems on closed manifolds. For instance, structurally stable dynamical systems on any closed 
surface have decompositions into invariant cells with the same dynamics, as we describe below in this 
paper; and their controllability follows from the same incidence relations. For bilinear systems these relations 
are extremely simple - they are given by linear inequalities, which is presented in this paper. 

More precisely, the main goal of this paper is to give sufficient conditions for controllability of bilinear 
systems satisfying the LARC (CCl condition) in three-dimensional case and have an open control set U (CC2 
condition). We do not consider yet the case n = 3 completely in all details, but restrict our attention to 
the "stable" or "general" systems (see [PDM]) as follows: we choose some finite subset U'^ of parameters 
Uk, k = I, N and assume (our stability or generality conditions CC3, CC4; sec below) that the systems of 
linear equations E^ for these parameters have different eigenvalues and their eigenvectors are different. 
Then we define - by a number of linear conditions on these eigenvectors - the open subsets CO'^{E^), which 
consist of geodesic triangles in the unit sphere. Finely, we prove that the projected system E^"" is controllable 
if some simple incidence relations are satisfied, see the Theorems A-C below. 



^For the sake of completeness note that all one-dimensional (n = 1) control systems x = A{u)x are controllable in our sense 
iflf A{u) ^ since all non-zero x have the same direction - the set of unit direction RP*^ is a one-point set. 
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Our arguments rely on the "Closed Orbit Controllability Criterium" from [M2], where it is proved that the 
control dynamical system on a compact manifold satisfying the LARC and with an open set of control vectors 
(see definitions in [M2]) is controllable if and only if through an arbitrary point goes some closed trajectory 
of the system. To make this paper self-content we reprove the two-dimensional version of this criterium here. 
Then we find sufficient controllability conditions such that compositions of shifts Pj: along trajectories of some 
finite number A{uk) provide us with the collection of closed trajectories going through an arbitrary q. 

The authors sincerely thank Professor D. Elliott for valuable comments. 

1. Bilinear control systems in 

As we mentioned above, the question if the given bilinear system E in R^ is controllable or not can be 
reduced to the same question about controllability of the projected system S^"" on the sphere S"^. In this 
section we consider first the more general question of the controllability of systems Em on closed (connected, 
compact without boundary) surfaces - we reprove for n = 2 the closed orbit controllability criterium from 
[M2]. Then, using this criterium, we give sufficient conditions for the controllability of the projected systems 
E^"' on the two-sphere S^. These conditions we formulate in terms of finite control subsystems of E^"' - for 
finite number of parameters Uk,k = 1, ...,N we consider systems 

(5) m = V{q{t),u,) (Ef), 

and find conditions on the set of E^'^ such that compositions of shifts along trajcctoricis of these systems 
(i.e., chronological products, see [M2]) give us closed trajectories going through all points q of the sphere. This 
is equivalent to say that every two points q and p may be connected by some trajectory of the system with 
piece-wise constant control functions with values Uk, k = 1, N. 

2. Controllability conditions 

The first condition we require is that the collection of vector fields Vk{q) = V{q,Uk),k = 1,...,N satisfies 
the LARC which for or any other surface is the following: 

CCl. In an arbitretry point q € the set of all control vectors {V{q)} contains two linearly 
independent vectors. 

For the control dynamical system on the surface the set V{q) = {V{q,u)\u € U} of all control vectors 
at the point q is the subset of the tangent plane TqM"^ to at the point q. The union V of all V{q) is the 
subset of the union of all planes TqM"^ - the tangent bundle TM"^. The second condition we require is 

CC2. The set V is an open subset in the tangent bundle TM"^. 

For projections E^"" of bilinear systems with an open control set U this second condition CC2 follows from 
the first one; i.e., from the LARC. 

Lemma 1. (Closed orbit controllability criterium for n = 2) Let Em be some control dynamical system on 

the closed connected surface IvP : 

(6) q{t) = V{.q{t),u). 

If the system Em satisfies CC1-CC2 above, then it is controllable if and only if through an arbitrary point q 
goes some closed and non-trivial orbit of the system Em • 

Proof. Let B be some coordinate neighborhood of the point q with local coordinate functions : B ~> R,i = 
1,2 such that x^{q) = 0. Then the tangent plane to at the point q' & B is generated by two coordinate 
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vectors ei{q'),i = 1,2 - evaluations at this point q' of the coordinate vector fields Cj = d/dxi, while for the 
control vectors holds: 

(7) V{q',u) = V\q',u)e,{q') + V\q',u)e2{q') 

for some functions V^{q', u),i = 1,2. If q{t),0 <t<Tis the closed trajectory of Em with the control function 
u{t) going through the point q, q = q{0) = q{T), then 

(8) ^^x\q{t))=V\q{t),uit)). 

Because CC2 is satisfied, for sufficiently small positive e all vectors V' from T^^^-^M for < t < e which are 
e-close to V{q{t),u{t)) belong to the control set V{q{t)). Therefore, for an arbitrary small vector v = (v^^v"^) 
with {v^Y + (^^)^ < solution of the equation 

(9) |4(i) = y^(g(i),«(i)) + ^^ 

also gives us the trajectory qv{t) = {x]j{t) , x1(t)) of the control system Sm- Since, for a given < t < e 
the set of all points qv{t) is an open set which contains q{t), we see that the set of points 0^{q) which are 
reachable from q contains some open neighborhood of (every) q{t) for < t < e. By standard theorems on 
ordinary differential equations, the image of every such open neighborhood under all shifts P'^ ,Q < t < T — t 
along trajectories of the system Em with the control function u{t) is an open-neighborhood of the trajectory 
q{t),Q < t < T. Therefore, the positive orbit 0~^{q) of the point q contains together with arbitrary closed 
trajectory going through q some its open neighborhood - which we denote by B^{q). The negative orbit 
0~{q) of q is the set of points from which we can reach the point q by trajectories of our system Em- It is by 
definition the positive orbit of the "minus"-EM control system, i.e., the system —Em on M with the control set 
—V{q). Easy to see that these two systems satisfy or do not satisfy our conditions CC1-CC2 simultaneously. 
Therefore, we may repeat our arguments and prove that 0~{q) contains with an arbitrary point q{t) for 
negative t < some its open neighborhood B~{q{t)). Because for our closed orbit it holds q(t) = q{T + t) 
we see that for arbitrary < r < T + i this neighborhood is an image of some open neighborhood of q{T) 
under the shift p^-^+* along trajectories of Em- Which implies that there exists some open neighborhood 
B^{q) of the closed orbit such that from its every point we can reach q. Then by definition the system Em is 
controllable on the intersection Bg = B~{q) D B^{q), which is again an open neighborhood of the closed orbit 
q{t),0< t < T. 

Now controllability of Em on the connected surface M follows from its compactness, see [M2]. For every 
q ^ M there exists some closed orbit going through this point. Every such orbit has an open neighborhood 
Bg on which the system Em is controllable. The union of all Bg covers M^. Since is compact there exists 
some finite sub-covering 

(10) M2 =Ui=l,...,c-Bg,. 

Because is connected for two arbitrary points q and p we may find a finite sequence Bg.^ , Bg.^ from the 
finite covering (10) such that q G Bg^^ , p G Sg.^ and Bg^^ PI -Bg.^^^ ^ 0. Take qo = q; some qs G Bg^^ n Sg.^^^ 
and qk+i = P- Since Em is controllable on each Bq^ there exists a trajectory from qs to qs+i for all s = 0, fc. 
Taking the composition of these trajectories we obtain the trajectory of Em going from q to p. Since q,p where 
arbitrary this proves that Em is controllable on M^. 

To complete the proof of our controllability criterium (... "and only if part) we note that if the system is 
controllable than for arbitrary q there exists the trajectory 7 of the system going from q to the point q{—e) 
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for some other non-trivial trajectory q{t), —e<t<eof the system Em going through q = q{0). Then the 
composition of 7 and q{t),—e < t <Q gives us the non-trivial closed trajectory through the point q. The 
Lemma is proved. 

Now we go back to the projected bilinear systems. 

For every fixed Uk the system S^'^ is the projection of the system of linear differential equations in E?: 

(11) x{t) = A{uk)x{t) (Sfc). 

Denote by \i{u) , \2{u) , X'^iu) eigenvalues of A{uk) and by Ei[u),E2{u),E^(u) the corresponding unit eigen- 
vectors. 

Here we will split our consideration into two cases: first is when all eigenvalues are real, and the second one 
- when one is real and other two are complex conjugate numbers Xq and Xq 

3. Real eigenvalues. Jordan cells decomposition and dynamics 

Here wc further assume 

CCS. The eigenvalues A^,i = 1,2,3 are different: 

(12) Ai < A^ < A^, 

which means that the linear operator A{uk) is stable - similar to all other operators which are close enough 
to it. Then in the basis E^(u),i = 1,2,3 the solutions of the system (11) arc given by 

(13) Xk{t) = {xlit) = xl{0)e^"\xUt) = 4(0)e^^,4(i) = 4(0)e^3*). 

We see that closer of domains in where coordinates have the same sign are invariant sets of this system. 
Correspondingly, triangles on the unit sphere with vertices iE^, i = 1, 2, 3 are invariant subsets of the projected 
system T,^^ . Sides of these triangles (also invariant subsets) are intervals of big circles (geodesies) in the unit 
sphere 5^. We denote by the circle around Ef consisting of such intervals; e.g., the circle C| consists of 
the intervals from E![ to then from i?| to —E^, then from —E^ to — i?| and from — i?| to . Inside each 
invariant triangle the position of the point q may be given with the help of coordinates x^- For definiteness 
we consider the triangle = A'^(-|- -|- -|-) with vertices Ei,E2,E^, where all x]. are non- negative. For an 
arbitrary point q inside this triangle its coordinates {x\{q) , x^q) , x'j^^q)) in the coordinate system Xk provide 
homogeneous coordinates {x^{q) : x1{q) : x]^{q)) for the unit direction q in the projective sphere RP^; it holds 
{xl{q)f + (xKq))^ + {xl{q)f = 1. The vector field Vk{q) = V{q,Uk) of the dynamical system S^'^ due to (12) 
is stable (sec [PDM]) and has six zeroes: two sources at iE'f , two saddles at ±i?2 and two sinks at ±£"1. 
Circles C| consist of saddle separatrix incoming to saddles ±£^|, and circles - of outgoing separatrix from 
saddles to sinks ±E^. 

If we denote by Pj: the shift along trajectories of T,^^ - the one-parameter group of diffeomorphisms of S"^ 
generated by the vector field Vk and given by solutions of (11) - then from (13) (considering mutual ratios of 
coordinates a;^ (t) as t — )• 00) we easily conclude the following description of the dynamic of such shifts inside 
our invariant triangle Afe. 

Lemma 2. 

1.1. For an arbitrary point q in the interior of and arbitrary small e-neighborhood Si(e) of E^ there 
exists a trajectory ofT,^ starting at some point q{€) in B^{e), going through q and converging to E^. 

1.2. For an arbitrary curve q{t),—e < t < e which is transversal to the outgoing separatrix £^2^3 '^^'^ 
intersect it at the point q = q{0), and another curve q{t),—e < t < e which is transversal and intersect the 
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incoming separatrix at the point q = q{Q), such that both q{t) and q{t) belong to the interior of the 

triangle Ak for positive t, i, for small enough positive t < e' there exists a trajectory of T,^^ starting at some 
point q{t) for some positive i = t{t) and going through the point q{t). 

Note that the claims of the Lemma 2 are valid also for non-stable dynamical systems when (12) is wrong, 
but equal eigenvalues correspond to A{uk) given by the Jordan form matrix on the subspace generated by 
corresponding eigenvectors. In other two cases we have or Vfc identically zero - when is diagonal and has 
equal eigenvalues; or Vk equals zero on some side of the triangle Afe corresponding to the two-dimensional 
subspace, where the action of A{uk) is diagonal with equal eigenvalues. 

4. Complete subsystems. Real eigenvalues 

Let = {uk\k = 1, AT} be some finite subset of the set U of control parameters, and T,^^ as above 
denotes the projection (11) of S^. We call complete if every two points q,p on the unit sphere can be 
connected by some trajectory of S^'' with the piece- wise control function u{t) : [0, T] — > U^; i.e., taking values 
in the finite subset U'^ . Now we give some sufficient conditions which guarantee that is complete. First 
we assume in addition to CCS that our subsystem is general: 

CC4. All eigenvectors are different, and E^ does not belong to any plane H^^ generated 
by Ei^ and Ei^ for all different si, 82,83 = 1,2,3. 

Next, for an arbitrary point q we generate two sets: first is the set CO^{q) of saddle points E2, and the 
second set CO^{q) - some open subset of the unit sphere, which is the union of (geodesic) triangles in this 
sphere with sides being big circles. Finely, our completeness criterium is that for some index s both Ef and 
—Ef belong to CO^(i?|) and CO^(— sec the Theorem A below. This criterium is given in terms of a 
number of simple inequalities, which describe mutual positions of systems of eigenvectors {Ef,i = 1,2,3}. 
The construction of the sets CO^ and CO^ proceeds as follows. 

Construction of CO^ and CO^ is by repetition of steps of two types. 
First type step. 

Take some point q in the interior of A;; and assume that some circle Cg "cut" this point from the vertex 
E^ of the triangle A^,. Then by the Lemma 2 above the trajectory q{t) = Pl{q) of T,^^ issuing from q converge 
to E3, and therefore, intersect the circle C3 at some point q = P^il)- Then ior t <i and close to t the points 
q{t) lie in some, say Aj(+ H — ) - the triangle with vertices ,£'2, —El^-, while for f > i and close to t the 
points q{t) lie in the triangle Aj(+ + +) - the triangle with vertices E{,E2,E'^. Applying the Lemma 2 again, 
now for the triangles Aj (+ H — ) and Aj(+ + +) we see that with the help of shifts along trajectories of TF^ 
we may approach arbitrarily close an arbitrary point of the semi-circle consisting of outgoing separatrix from 
E2 - the union of intervals —E^^E2 U E^^E^^, which we denote Sep2{+)- Since our trajectories under the flow 
Pj are coming from a neighborhood of the point q, they converge to the separatrix only from one side, the 
points which are reachable from q{t) always belong to the same half-sphere relative to Cg as the point q. Note 
in addition, that for the given point q in the interior of A^ and separated from E^ by Cg we may find some 
open neighborhood B{q) such that the same is true for all q' € -B((?). Thus, from the Lemma 2 we conclude 
the following. 

Lemma 3. If the point q in the interior of Ak is separated from E^ by C|, then for all q' in some small 
neighborhood B{q) of q, and an arbitrary point p on Sep2(+) which is not end point of it, i.e., p 7^ i-E'a there 
exist an open neighborhood B{p) such that all points in B{p) are reachable from q' if they belong to the same 
half-sphere relative to Cg as the point q. 
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Let us write down our conditions in an analytic form. If S{+) and >S'(— ) denote the two half-spheres in 
which the big circle C divides the unit sphere, then the point q belongs to one of them depending on the sign 
of the scalar product of q with the vector normal to the plane in which this circle lies. Or, if C lies in the 
plane n gcncirated by Ei and E2 - the vector normal to 11 is the vector product [Ei , E2] , and the points q and 
p belong to different half-spheres if and only if the corresponding products have different signs: 

{q, [EuE2])ip,[EuE2])<0, 

or in a more symmetric form: 

det{q, Ei,E2)det{p, ^1, ^2) < 0, 

where by det{X, Y, Z) we denote the mixed product of three vectors X, Y, Z. Define the function: s{q, C,p) = 1 
if (14) is satisfied, and s{q, C,p) = —1 if not. Then the points q and E^ are separated by C3 if 

(14) s{q,Ci,E'^) = l. 

The triangle is the intersection of three half-spheres out of six S'l{±), i = 1, 2, 3 in which the big circles 
C|^,i = 1,2,3 - its sides - divide the unit sphere. So, the condition that the point q belong to the interior of 
Afc is that the vector q belongs to the same half-sphere relative to as the vertex E^ (i.e., which does not 
belong to this circle): 

{q, [El El]){El [El E^]) > 0, (5, [EIE^]){EI [E^ E^]) > 0, {q, [E^ EI){EI [E^ E^) > 0, 

or 
(15) 

det{q, ElEl)detiEl,El E^) < 0, det{q, ElE^)det{El E^ eI) < 0, detiq, E^ E^detiEl E^ e1) < 0. 

Again, we define the function d{q, A) = 1 if the point q belongs to the interior of the triangle A = A{Ei,E2, E3) 
with vertices Ei,E2, E3; and d{q, A) = — 1 - if not. From the above we see that d{q, A{Ei,E2, E3)) = 1 if and 
only if 

(16) s{q,Q,E,)^-l for all i = l,2,3, 

where Ci denotes the big circle in the sphere going through the points Ej,j 7^ i. 

If the incidence condition (15) holds for some k, and for this k the separation condition (14) holds for some 
other j, then we set: 

(17) EiGCO'iq). 
Second type step. 

Let the point E2 S CO^{q) belongs to some triangle A;, say Ai(-|--|--|-), such that the semi-circle Sep2{+) cut 
from this triangle some vertex El; i.e., such that end points ±£^3 of 5ep2(+) lie outside A;. Then, depending 
on s = 1,2,3 the shifts Pf along trajectories of of the points p from the chord Tj = Sep2{+) H A; of 
the semi-circle 5*6^2 (+), which is inside the triangle A(, cover some open subregion ilf, which we include in 
CO'^{q) (and define it precisely in a moment below). Since such points p belong to the interior of A; every 
trajectory Pf (p), < f ' < T" can be continued over some open interval — e" < i" < such that the extended 
trajectory contains points from both sides of the chord . Because all points p from F^ are interior points 
of the separatrix 5*6^2 (+) (end points ±El of 5*6^2 (+) lie outside A;) we conclude from the Lemma 3 above 
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that all points of the open domain are reachable by trajectories from q and from all points q' of some its 
small neighborhood B{q). 

In addition, when Fj cut the vertex E^^ for s 7^ 3 - it has some end point on the incoming separatrix E\E2. 
Therefore, arguing as in the Lemma 3, we also include the limit point of this separatrix E2 in CO^ under this 
step of the second type. 

Write down the analytic form of our arguments. The point E2 belong to Ai{E[E2E^^) if 

(18) d{Ei,Ai{E[EiEl)) = l. 

The endpoints of the semi-circle Sep2{+) do not belong to this triangle if 

(19) d{±Ei,Ai{E[E'2Ei)) = -l. 
Then the chord cut from this triangle the vertex E[ if 

(20) s{E[,Ci,Ei) = l and s{E[,Ci, Ei,) = 1. 

In this case the union of all trajectories of Sj"^ issuing from points of F^ covers the part of A; lying in the 
same half-sphere relative to C3 as E^: 

(21) nj = {p\d{p, Ai) = 1, sip, Ci,El) = -1}. 
If the chord F; cut from this triangle the vertex Eji i-^-) 

(22) s{ElClEi) = l and s{ElClEi) = l, 

then the union of all trajectories of issuing from points of F^ contain the part of A; lying in the same 
half-sphere relative to Cg as E2 - small geodesic triangle with the vertex E2 and the opposite side Fj . Indeed, 
an arbitrary point p inside this triangle by the Lemma 2 lies on some trajectory p{t) coming from some small 
neighborhood of E{ - outside this triangle, and therefore intersecting Fj at some point. Hence, the points of 

(23) = {p\d{p, Ai) = 1, s{p, Ci,E'2) = -1}. 

are reachable from F^ at this case (22). If the chord F^ cut from this triangle the vertex E^; i.e., 

(24) siE[,Ci,Ei,) = l and siE^Ci, E',) = 1, 

then similarly all the points from the small triangle with the vertex E^ and the opposite to it side F; are 
reachable from the points of this side, and 

(25) ni = {p\dip, AO = 1, s{p, ClEi) = -1}. 



After this step we add to CO'^{q) the open region Q,j and to the list CO^{q) the vertex E2 if (20) or (22) 
holds. The result of the step of the second type we formulate as follows. 



DIRECTION CONTROL OF BILINEAR SYSTEMS. I. 



9 



Lemma 4. Let E2 G CO^{q) belongs to some triangle A;, and the semi-circle separatrix Sep2{-\-) cut from 
this triangle some vertex E^.; i.e., such that end points ±i?| of Sep2{+) He outside A;. Then, depending 
on s = 1,2,3 the open set flj defined by (21) if s = 1, (23) if s = 2 or (25) if s = 3 is contained in the 
positive orbit 0^{q') for all points q' from some small neighborhood B{q) of q. If s = 1,2 then in addition the 
conclusion of the Lemma 3 is true for all interior points p of the separatrix Sep2{+)- 

Now we give our sufficient controllability criterium for the subsystem of T,^^ , k = 1,...,N where all eigen- 
values are real. 

Theorem A. Let Y^p^ be some projected system which satisfies the LARC ( CCl ), the openness condition 
CC2) and generality conditions CC3, CC4- Assume that for some finite subset = {uk\k = 1,...,N} of U 
and some index s after some finite number of steps, each of which is of the first or second type as above, it 
holds 

(26) El and - E{ belong to CO'^{E^) and CO'^{-EI). 

Then the system S?"' is controllable. 

Proof Take an arbitrary point q. As in the lemma 1 we see that, from CCl and CC2 it follows that both the 
positive and negative orbits of q contain some open sets - neighborhoods B^{q{e)) and B~ {q{—e)) of some 
trajectory q{t) going through q. Thus, taking partition of the unit sphere into the set of triangles As(±±±) we 
sec that there is some of these triangles, denote it by A+ with the vertex i?| or — i?| - say £^1 for dcfiniteness; 
which intersect the positive orbit of q, say at some point q' . Therefore, by the lemma 2 there exists some 
trajectory from q' to some p' in an arbitrary small neighborhood of this vertex E^. If we take this small 
neighborhood to be B{E^) as in the Lemmas 3 and 4 above, then wc find a trajectory from q to some point p' 
from which we can reach all the points of CO^{E^) and and —£{ in particular. Since CO^{E^) is an open 
set, we can reach all the points from some small neighborhoods of Ef and —Ef, which we denote B'{Ef) and 
B'i-Et). 

Since the negative orbit of q also contains some open subset, it has non-empty intersection, say the point 
q", with the interior of some of the triangles A J with the vertex Ef or —Ef, say Ef for dcfiniteness. Again, 
by the Lemma 2 the point q" belongs to some trajectory issuing from some point p" which we may choose in 
an arbitrary small neighborhood of this vertex. Finely, the combination of these orbits: from q to q' , then to 
p' , then to p", then to g" and back to q gives us the closed orbit going through q. Since we choose the point 
q in an arbitrary way, we have proved that the control dynamical system S*^, /c = 1,...,N satisfies the closed 
orbit controllability criterium and is controllable by the Lemma 1. 

The Theorem A is proved. 



5. Complete subsystems. Complex eigenvalues 

When the system Ej has one real eigenvalue \\ = A'^ and the pair of complex-conjugate eigenvalues 
A2 = Re{\-'(j) + y/^Im{X(^) and A3 = Re{XQ) — ■\/^Jm(A;^) then its solutions are given by Xj{t) = 
{x](t),x]{t),x^j(t)) where 

x]{t)=x]{0)e^^' 

and 

x^j{t) = x'^j{0)cos{Im{X{.)t)+x^j{0)sin{Im{X{.)t))e^<^c)t, 
:r3(<) = -x^.{0)sin{Im{X^)t) + x'^j{0)cos{Im{X^)t))e^^^^o)t. 
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Here we have only two two-dimensional invariant subsets - half-spaces = {xj > 0} and R^ = {x^ < 0}. 
Correspondingly, the projected system T,^^ has two invariant cells which are half-spheres. They intersect along 
the circle which is also an invariant set of S^' . Another invariant set - two zeros ±El of Vj. If 

(27) A], > Re{X^), 

then zeroes ±El are sources of the dynamical system E^*^, while is an attracting cycle, and if 

(28) A], < Re{X^), 

then zeroes ztE^ are sinks of the dynamical system S^"^, while is a repulsive cycle. The analog of the 
Lemma 2 describing the dynamic of T,^^ is the following. 

Lemma 5. 

1. If El is the source, and is an attracting cycle, then for an arbitrary point q in the interior of Aj and 
arbitrary small e-neighborhood B{El,e) of Ef there exists a trajectory q(t) ofYF^ starting at some point q{e) 

in B{E{,e), going through q and converging to the big circle : for an arbitrary point p of this circle and 
arbitrary q there exists a sequence tn{q,p) — >■ oo such that q{t„) converge to the point p. 

2. Let E\ be the sink, - the repulsive cycle, and p{t) some curve intersecting transversally at the 
point p = p(0). Then for an arbitrary point q in the interior of Aj and arbitrary small e-neighborhood B{p, e) 
of p there exists a trajectory qi{t) ofEY starting at the point p(t) in B{p,e), going through q and converging 
to the sink El . 

Using arguments from the Theorem A above we see that the LARC (CCl) and the openness (CC2) 
conditions with the help of the Lemma 5 immediately imply that the positive orbit of the source El (under 
(27)) contains some open neighborhood of the half-sphere A-' (-l-) (and 0^{—El) D A-' (— ) correspondingly); 
while if El is a sink (when (28) holds) even more is true: the positive orbit of any small curve p{t) intersecting 
transversally at this point p contains the whole unit sphere. 

Having at least one such E^' in our finite subsystem simplifies our controllability considerations: Let 
denotes the parameters Uh in U'^ such that Sfe has only real eigenvalues, and CO^{q), CO'^{q) constructed 
as above for the subsystem . Then repetition of the arguments in the proof of the Theorem A give us the 

following. 

Theorem B. Let "S'"' be some projected system which satisfies the LARC ( CCl ), the openness condition CC2) 
and generality conditions CCS, CC4; andU^ = {uk\k = 1,...,A'^} some finite set of parameters such that for 
one of them j the system has complex eigenvalues. Assume that for the subset = {uk\k = 1,...,A'^} of 
parameters such that all eigenvalues of Sfc are real, and some index s after some finite number of steps, each 
of which is of the first or second type as above, it holds: 

(29) El and - El belong to CO^{El) and CO'^{-El), 
if (27) is satisfied; or 

(30) 5*^ intersects both CO^{EI) and CO'^{-EI), 
when (28) holds. Then the system S^'' is controllable. 

The most simple case is when we have two systems Ej and Sigmai with complex eigenvalues with different 
types of dynamics: when has the attracting cycle , while E^' has the repulsive cycle 5'. Then, using 
the system S^'" with attracting cycle we may construct a trajectory from an arbitrary point q to an arbitrary 
point of the attractive cycle - say, to its intersection with the repulsive cycle S*', and then continue this 
trajectory further to an arbitrary point with the help of the system E^*^. We have the following controllability 
condition. 
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Theorem C. Let T,p^ be some projected system which satisfies the LARC (CCl), the openness condition 
CC2) and generality conditions CCS, CC4; and in the finite set = {uk\k = 1,...,N} of parameters for 
one of them Uj the system has complex eigenvalues satisfying ( 27), and for another Ui the system has 
complex eigenvalues satisfying (28). Then the system is controllable. 
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